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Abstract 

The leading m# radiative correction to the Higgs decay width into 
a pair of weak vector bosons is calculated at the two-loop level, using 
the equivalence theorem in Landau gauge. The result indicates the 
breakdown of perturbation theory if the Higgs boson is heavier than ~ 
930 GeV, in spite of the smallness of the one-loop radiative correction. 



1 Introduction 



After the recent discovery of the top quark at the Tevatron, the mechanism 
of spontaneous electroweak symmetry breaking undoubtedly remains the 
most obscure part of the standard model. The screening theorem |l], 
which states that the leading contributions in m# cancel in the radiative 
corrections to low energy processes, results in the present loose limits on the 
Higgs mass derived from radiative corrections. 

A Higgs mass of the order of 1 TeV is definitely not excluded by the 
present data. A recent analysis of the LEP results even points in the direc- 
tion of a heavy Higgs ||]. However, other analyses give somewhat different 
results The uncertainty is apparently too large to reach a conclusive 
estimate of the Higgs mass. 

A large Higgs mass means strong interactions in the spontaneous sym- 
metry breaking sector, and would challenge the theory because of the break- 
down of the perturbative approach. For a not too heavy Higgs, one may try 
to resum the asymptotic perturbative series by means of Pade approximants 
or by using various nonlinear sequence transformations such as the Levin's 
sequence transformation. Such techniques proved useful on simple quantum 
mechanical examples Q, but strictly speaking there is no proof that they 
would converge in the case of radiative corrections in the standard model. 

At the same time, strong couplings in the spontaneous symmetry break- 
ing sector may result in a rich, nonperturbative spectrum of new phenomena 
at the TeV energy scale. In the limit of infinite Higgs mass, the Higgs sector 
reduces formally to the Lagrangian of the nonlinear sigma model, with the 
role of the chiral symmetry played by the custodial symmetry. In the low 
energy limit one can use chiral perturbation theory to derive anomalous self- 
interactions of the electroweak gauge bosons ||. At the TeV energy scale, a 
nonperturbative spectrum of resonances may appear by analogy with pion 
physics [j?]], and restore partial wave unitarity in longitudinal vector boson 
scattering. The BESS model was originally proposed as a possible descrip- 
tion of such bound states of the electroweak symmetry breaking sector ^ . 

One can estimate the validity range of perturbation theory from unitar- 
ity violation of tree level amplitudes || [l(J, and from the size of one-loop 
radiative corrections [111], 12]. A perhaps more explicit criterion to estimate 



the range beyond which the asymptotic series stops converging at all re- 
quests knowledge of radiative corrections at two-loop order. The value of 
each term in the perturbative expansion of an amplitude depends on the 
renormalization scheme (see for instance ref. [13]), and therefore the range 
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of validity of the perturbative approach is scheme dependent. However, for 
a given scheme, at the point where the two-loop corrections are as large as 
the one-loop ones, the perturbative series shows no sign of convergence, and 
there is no warranty that the tree or the one-loop level amplitudes calculated 
in that scheme are reasonable approximations. 

Calculating two-loop radiative corrections in the electroweak theory is 
not an easy task. However, the problem becomes considerably simpler if the 
external momenta can be neglected. It was shown by van der Bij and Velt- 
man that any two-loop Feynman diagram evaluated at vanishing external 
momenta can be expressed in terms of the derivatives of the so-called mas- 
ter diagram, for which they found an analytical formula in terms of Spence 
functions [14]. This made it possible to evaluate two-loop heavy Higgs ef- 
fects on low energy observables such as the p parameter the masses of 
the vector bosons and the trilinear selfcouplings of the vector bosons 
[16|. As expected, these effects are rather small because of the screening 
of the heavy sector, and perturbation theory only breaks down for a Higgs 
particle as heavy as 3 — 4 TeV. 

To calculate two-loop corrections to processes involving the symmetry 
breaking scalars, one has to evaluate two-loop diagrams with finite exter- 
nal momenta. This is a more difficult problem, and calculations of physical 
quantities have not been performed until recently, in spite of the consider- 
able effort which was devoted to solving massive two-loop integrals. The 
reason is that in general the massive two-loop scalar diagrams with arbi- 
trary masses and external momenta cannot be expressed in terms of known 
and easy to evaluate functions, like polylogarithms. For instance, the two- 
loop selfenergy diagram with three propagators was shown to be related to 
a certain generalization of the hypergeometric series, the Lauricella function 

A number of techniques were proposed to deal with some two-loop dia- 
grams. Without exhausting the list, let us merely mention a few. Analytical 
results exist for certain diagrams evaluated at special values of the masses 
and of the external momenta Jig, Hjj]. Several integral representations were 



proposed for some diagrams po[ |22| , [2l| , 23], as well as Monte-Carlo inte- 
gration over the Feynman parameters |p8| , |29f| . Momentum expansions were 
worked out too ]^4|, ^q] , and conformal transformations combined with 
Pade approximants or Levin's sequence transformation were used to extend 
their convergence domain [27]. 

This paper uses a recently proposed method to deal with massive two- 
loop diagrams in a systematic way [30]. It can be used, at least in principle, 
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to evaluate any massive two-loop scalar integral with controllable accuracy, 
and is suitable for implementation in a computer program to generate and 
calculate automatically the Feynman graphs relevant for a given physical 
process. The same method was used to calculate the Higgs selfenergy at two- 
loop level and to extract the leading corrections to the location of the Higgs 



resonance as seen in fermion scattering [30]. This gives a perturbative bound 
on the Higgs mass of ~ 1.2 TeV. The two-loop heavy Higgs corrections to 
the Higgs fermionic width give a somewhat lower perturbative bound of 
~ 1.1 TeV ||, H§. 



This paper presents the two-loop corrections to the H — ► Z°Z° and H — ► 
W + W~ decays in the heavy Higgs limit. More precisely, only the leading 
ma corrections will be retained, which grow like mjj- at two-loop level. As 
expected, these processes allow one to derive a stronger perturbative bound 
on the mass of the Higgs boson. The two-loop correction becomes as large 
as the one-loop one for ran ~ 930 GeV, indicating that, at least in the OMS 
renormalization scheme, the perturbation theory is not reliable anymore. 

2 Lagrangian and renormalization 

Since one is only interested in radiative corrections at leading order in the 
Higgs mass, the natural choice is to use the equivalence theorem and the 
Landau gauge. 

The equivalence theorem [32[ — [^] relates the Green functions with ex- 



ternal longitudinal vector bosons Vl to the Green functions with the corre- 
sponding Goldstone bosons (f> replacing the vector boson legs: 

A(V£ ,V[ 2 ,...,vt) = {tC) n A{^ ,<?*,..., tfh) + 0{^r) . (1) 

Such relations are a consequence of the Slavnov-Taylor identities of the 
theory. The coefficient C with which the amplitude has to be multiplied for 
each external longitudinal vector boson replaced by a Goldstone boson is 
gauge dependent and is in general not unity beyond tree level 35]. 



The calculation is simpler if performed in Landau gauge. In this gauge 
the coefficient C in eq. 1 is one at leading order in ran fl35|| , and the dia- 
grammatics is considerably simplified. There are no Goldstone boson-vector 
boson mixing terms, and diagrams containing fermions, vector bosons or 
Fadeev-Popov ghosts need not be taken into account since they are non- 
leading. The only diagrams which give leading contributions in ran are 
those which contain only the symmetry breaking scalars. 
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The scalar sector of the standard model reads: 




where Hq, an d z o denote the bare Higgs and Goldstone fields. The 
gauge coupling constant g can be defined by using the muon decay as g 2 = 
Ay/lmly Gf, with G F = 1.16637 • 1CT 5 GeV~ 2 , and m w = 80.22 GeV. 
The energy scale where the gauge coupling constant is defined is actually 
irrelevant within this context because g is not renormalized at leading order 
in rriu- St is the tadpole counterterm which will be fixed by the condition 
that the vacuum expectation value of the Higgs field, v , should not receive 
quantum corrections, and is related to the selfenergy of the Goldstone bosons 
at zero momentum transfer []36| . 

To renormalize the theory, one splits the bare Lagrangian of eq. 2 into 
the renormalized Lagrangian and counterterms: 



We use an on-shell renormalization scheme with field renormalization. 
In this renormalization scheme, the counterterms can be determined from 
the following renormalization conditions: 




W 5 



(3) 



^HH(k 2 = m 2 H ) + i 5m 2 H — i5t + i 5m 2 H 5Zh — iSt 5Zh = 

d - 

-^Y< HH {k 2 = m 2 H ) + i 5Z H = 
^w+w- (k 2 = 0)-iSt-i St 5Z G = 



—t w+w -(k 2 =0) + iSZ G = 
T, w+W - (k 2 = 0) + i Sm^r = , 



(4) 
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where E contain the loop and loop-counterterm selfenergy diagrams, but 
not the pure counterterm diagrams. T iw + w -(k 2 ) is the coefficient of the 
—g^u piece of the vector boson selfenergy. Because of the Ward identity 
dm 2 ^ = —m^ySZc , it's not really necessary to consider the gauge sector 
to calculate the mass counterterm of the vector bosons, as done in the last 
line of eqns. 4, but this was done nevertheless because it provides a useful 
check on the calculation. Another check which was done is to compute 
the Higgs tadpole diagrams and to check that they cancel upon the tadpole 
counterterm derived from the Goldstone boson selfenergy at zero momentum 
transfer. 

At one-loop level, one needs to evaluate the H-H, w-w and W-W selfen- 
ergies whose topologies are shown in fig. 1, to find the following counterterms 
at order q 2 



xJl-loop) _ n 2 ™H_ /J™H_W2 j™H_ r_ jj_ 3 _ 3_7 

° ~ 9 - 9 Wit* 1 4e + 8 8 
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+_ L6 8 + 2 8 " 

r 7 (i-ioo P ) 2m 2 H( m 2 H /2 1 1 /3 7 \, 

^ G " 5 <: ( 4V } l6^ i_ 8 +e V 32"l6j j • (5) 

In these expressions the space-time dimension is n = 4 + e, and C7 
is the Clausen function, Cl(x) = Y^=i s i n {nx)/n 2 . Numerically, C7(|) = 
1.0149416064 .... The one-loop counterterms are needed at order e because 
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these terms combine with the poles of one-loop scalar integrals to give finite 
contributions at two-loop order. 

By evaluating the two-loop H-H, w-w, and W-W selfenergies, one finds 



2 

the following 0{{g 2 — #-) 2 ) counterterms pO, 31 
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4vr / u 2 7 (16tt 2 ) 5 

+-^-(169 - 967 - 24^71-) 

-(4.785031 ± 4.2 • 10~ 5 ) - i (0.412438 ± 1.6 • 10~ 5 ) ] } 
(2-ioop) _ , 2 m 2 H 2 m 2 1 3 1 13 

^ G " 15 m 2 y M 4vr / . 2j (16vr 2 ) 2 [ 32 e 128 + 32 7 

192 64 16 V 3 ;J 

, y (2-loop) _ p r / 2 m H ^2 / m g 1 r 3 1 

-(0.62296 ± 2.5 • 10~ 4 ) - i (1.00233 ± 2.5 • 10" 4 ) ] } , (6) 

These counterterms agree with an independent calculation of the Higgs 
and Goldstone selfenergies by P.N. Maher, L. Durand, and K. Riesselmann 



3 The calculation 

The partial decay widths of the Higgs boson into a pair of vector bosons are 
given at tree level by: 
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(7) 



The one-loop radiative corrections at leading order in m,H can be calcu- 
lated from the diagrams shown in fig. 3. The resulting Hw + w~ coupling, 
including the order g 2 m 2 H /m 2 v corrections, is: 



19 5vr 2 3\Z3tt 
16 + ~48 8 

(8) 

and the partial widths of eq. 7 correspondingly get a correction factor 

167r 2 m 2 w { 8 24 4 J ' { ' 

The real part of the correction of eq. 8 and the corresponding correction 
to the widths of eq. 9 agree with the results of ref. |i~2f| . Note that at 
two-loop level one needs the imaginary part of the one-loop radiative cor- 
rection of eq. 8 as well because it gives a correction to the widths of order 
(g 2 m 2 H /m 2 w ) 2 . 

One problem related to the use of the Landau gauge is the presence 
of massless particles in the theory. Some Feynman diagrams may display 
mass singularities, and there is the problem of the arbitrariness of integrals 
of the type / d n p (1/f? 4 ) in the framework of dimensional regularization. 
Such problems do not appear in this calculation at one-loop level, since 
the only place where one encounters a / d n p (1/p 4 ) type integral is the 
kpkv piece of the vector boson propagator. They are however present at 
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two-loop level, and one needs a regularization procedure to deal with these 
problems in a consistent way. The approach adopted in this paper is to 
work in a nearly-Landau gauge, that is, to keep a small gauge parameter £ 
during the calculation, and to let £ — > in the final results. This amounts 
to giving the Goldstone bosons a small mass Tnyy. This procedure is 
consistent with discarding the diagrams containing fermions, gauge bosons 
and Fadeev-Popov ghosts on the internal lines, since they do not give rise 
to contributions of 0((g 2 w-^/^iy) 2 ) in the limit £ — ► 0. In this way, an 
additional check on the calculation is present, namely the cancellation of 
the poles and logarithms of the gauge parameter in the final result. This 
cancellation involves both the analytical and the numerical parts of the 
calculation. Another possibility, advocated in ref. pS(| , is to calculate the 
singular diagrams with vanishing Goldstone boson masses but at off-shell 
momentum, and to check that their sum remains finite when the external 
momentum is put on-shell. 

Let us now turn to the actual two-loop calculation. The main task is 
to calculate the two-loop proper diagrams corresponding to the Hw + w~ 
vertex. The topologies which are involved are shown in fig. 4. 

Considering the large number of Feynman diagrams which need to be 
calculated and the lengthy expressions which are involved at intermediary 
steps, the whole calculation was done by computer. The algorithm is essen- 
tially the same which was used to derive the two-loop counterterms given 
in the previous section. 

As a first step, a computer program generates all relevant Feynman di- 
agrams. This was done by giving by hand the possible topologies of the 
proper vertex diagrams, which are shown in fig. 4. For each given topology, 
the program then substitutes for each internal line all possible propagators, 
that is, H-H, z-z, w + -w~ , and w~-w + . This way many diagrams are gen- 
erated which contain nonexistent vertices. The program then compares the 
vertices of the diagrams generated with a complete list of the vertices of the 
theory, discards the spurious diagrams, and substitutes the actual expres- 
sions of the vertices in the correct diagrams. The combinatorial factors are 
automatically correct, provided one divides each topology by its symmetry 
factor. If a certain topology has m sets of equivalent lines or nodes which 
contain ni,ri2, • • • , n m elements, then its symmetry factor is n\\ n?}. . . . n m \. 

In the next step some algebra is necessary to bring the resulting Feynman 
diagrams into a standard form |3(J. After doing all the partial fractioning 
which is possible, the program decides which diagrams can be calculated ana- 
lytically and which ones need numerical integration. The two-loop diagrams 
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evaluated at vanishing external momentum can be expressed analytically in 
terms of Spence functions, and in this case their mass expansion is needed 
in order to retain only the 0((g 2 m? H / 'm^) 2 ) terms in the final expression. 
In the case of the two-loop diagrams evaluated at finite external momenta, 
the propagators with the same loop momentum are combined by means of 
Feynman parameters, and the resulting scalar integrals are calculated in 
terms of two basic functions, g and /, whose definitions are given in the 
Appendix. 

Some care is needed when introducing Feynman parameters in trian- 
gular diagrams. It is useful to parametrize in the same way the diagrams 
with a similar structure in order to avoid the unnecessary proliferation of 
expressions resulting from introducing Feynman parameters. There are two 
types of triangular diagrams which appear in this calculation, and they can 
be parametrized in the following way: 



1 _ T(ai + a 2 + a s ) 

[(p -h) 2 - M 2 ]«i [(p + k 2 ) 2 - M 2 ]«2 [p 2 - m 2 ]«3 ~ r(ai)r(a 2 )r(a 3 ) 

x r ix r iy 'W-.-Mi (10) 

Jo Jo [(p + kx) 2 - m\ + i 7? ]«i+«2-i-«3 

1 _ T(a 1 +a 2 +a 3 ) 

[{p - h) 2 - m 2 ]«i [p 2 - M 2 ]«2 [(p + k 2 ) 2 - m 2 ] a 3 ~ T( ai )T(a 2 )T(a 3 ) 



x 



f 1 dx f 1 d xjxyr^lxjl-y^^il-x)^- 1 
Jo Jo [(p + k 2 ) 2 - m 2 , + ir]] a i-+ a 2 +a 3 



where k\ and k 2 are the momenta of the external Goldstone bosons, M = 
mu and m = y/£ m\y are the masses of the Higgs and Goldstone bosons, 
and 



k\ = —xyki+x(l — y)k 2 

k 2 = -xy ki + (1 - x) k 2 

m\ = x[l - xy(l - y)\ M 2 + (1 - x)m 2 

ml = x [l-y(2-x)]M 2 + [I - x{\ - y)]m 2 . (12) 

The diagrams containing the second structure have at least one two- 
particle cut, and correspondingly an imaginary part. After introducing 
Feynman parameters, this translates into the presence of poles which one has 
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to avoid by a suitable choice of the integration path. With respect to the nu- 
merical integration over the Feynman parameters x and y, the parametriza- 
tion given in eq. 11 has the advantage that the solution of the equation 
m2 = has a simple structure in the integration domain (0, 1) x (0, 1), and 
therefore it is easy to find an integration path to avoid this singularity. 

Also related to this parametrization is the problem of using the Landau 
gauge to calculate the diagrams of fig. 5. Such diagrams taken separately are 
not well defined for physical momenta because of an endpoint singularity, 
and it is only their sum which is finite. More precisely, these diagrams 
have two two-particle cuts, and after introducing two Feynman parameters, 
one obtains an expression of the type in eq. 11. The singularities of the 
integrand in the xy plane lie on a curve given by 777,2 = 0, which intersects 
the boundaries of the integration domain (0, 1) x (0, 1). These singularities 
can be avoided easily by deforming the integration paths in a convenient way 
inside the square (0,1) x (0,1), but not on its frontiers, and therefore the 
integral is logarithmically divergent. One way to regularize this divergency 
is to calculate these diagrams slightly off-shell, for instance by keeping a 
small but finite 77 in eq. 11, and to take the limit 77 — > in the sum of these 
three diagrams. The problem is that also the Landau gauge limit must be 
taken at the end of the calculation, and this must be handled with care 
because these two limits do not commute. This is due to the presence of 
pieces which are nonleading in 777/7/ and ought to vanish in Landau gauge, 
but which display also the endpoint logarithmic singularity. Therefore, the 
correct way to take the limit is first to set £ — > 0, and then to go on-shell. 

A simple way to get rid of this problem is to introduce explicitly in the 
diagram of fig. 5 c) the terms £ rriyy 5Zq w + w~ and 1/2 £ 5Zq zz, which 
in fact exist in the Lagrangian of eq. 2 when one is not exactly in Landau 
gauge, but which only generate contributions which vanish when £ — > 0. 
Taking these terms into account in the diagram 5 c) ensures the cancella- 
tion of the nonleading contributions which have the endpoint singularity, 
after which the order in which one takes the limits £ — > and 77 — > be- 
comes irrelevant. Since the endpoint divergency is only logarithmic, no large 
numerical cancellations appear. Moreover, the cancellations among the di- 
agrams already occur at the level of the integrands, before the numerical 
integration is carried out, because the diagrams were parametrized in the 
same way. The result is thus numerically stable. 

Some checks on the algebraic part were also included in the program. 
Where possible, the analytical cancellation of the poles and logarithms of 
the gauge parameter was checked. Some subsets of diagrams must be free 
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of ultraviolet divergencies, in agreement with Bogoliubov's proof of renor- 
malizability. This is the case with the combinations c + d + r, k + l + s, and 
g + h + i + j + o, where the notations of topologies are defined according to 
fig. 4. Note that in order to check this, one needs first to set the one-loop 
field renormalization counterterms to zero. 

The algorithm described was encoded in a FORM [37] program. It takes 
approximately seven hours to generate all Feynman diagrams, and to per- 
form the necessary algebra and the checks on a NeXT computer. In the 
end, one obtains an analytical part and a number of numerical integrals 
over Feynman parameters of g functions, which were encoded in FORTRAN. 
Their evaluation took approximately 10 hours on an Apollo 9000/720 work- 
station. 

In the end, one obtains the following result for the Hw + w~ coupling, 
including the order (g 2 mfj/rriyy) 2 radiative corrections: 



■ 9 m 2 H 



2 mw 



1 + 



g 2 m 2 H 
16 7r 2 mly 



19 5vr 2 3\Z3tt 

H — — 7, \~ i7r 







( 92 


m%\ 2 


\IQ-K 2 


m w) 



16 48 
(.53673 ±4.1 • 10" 4 



log 2 



-i (.32811 ±3.1 • 10~ 4 )] } . 



(13) 



Correspondingly, the partial widths of eq. 7 get the following correction 
factor: 




This correction factor is shown in fig. 6 as a function of run- The 
two-loop correction has the same sign as the one-loop correction, and for 
mn ~ 930 GeV it becomes as large as the latter. This is an indication for 
the validity range of perturbation theory. This is a rather surprising result, 
taking into account that the one-loop radiative correction is quite small for 
such a Higgs mass, at the level of ~ 13%. 
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m H 
[GeV] 


T(H -» W+W~) 
[GeV] 


T(H -» Z + Z~) 
[GeV] 


r(H -> «) 

[GeV] 


total 
[GeV] 


400 


16.97 


7.920 


2.146 


27.04 


500 


36.82 


17.60 


10.89 


65.32 


600 


68.43 


33.16 


20.15 


121.7 


700 


115.7 


56.55 


29.09 


201.4 


800 


184.0 


90.38 


37.59 


311.9 


900 


280.5 


138.3 


45.58 


464.4 


1000 


415.4 


205.4 


52.92 


673.7 



Table 1: The main decay channels of a heavy Higgs boson at two-loop order 
(m t = 180 GeV). 

The partial decay widths corresponding to the main decay channels of 
the Higgs boson, including the two-loop 0((g 2 ™ 2 h / m w) 2 ) radiative correc- 
tions, are given in table 1. The tt channel, which was calculated with similar 
methods in ref. [ )31| j , is also given. It should be noted that by multiplying 
the tree level decay rates in eqns. 7 by the correction factor of eq. 14 some 
subleading terms are also generated. They start with 0(g 2 ) terms in the 
one-loop correction, and with 0(g 4 T r n? H /m 2 v ) terms at two-loop. Such con- 
tributions are of course incomplete, but they were not explicitly subtracted 
from the numerical results given in table 1 because they are formally of the 
same order as the theoretical uncertainty due to the full subleading contribu- 
tions they are part of, and also numerically negligible. A similar discussion 
holds for the tt channel as well. 

4 Conclusions 

The decays H — ► W + W~ and H — > Z°Z° were calculated at two-loop level 
in the limit of large Higgs mass. The calculation was performed in Landau 
gauge and by using the equivalence theorem, in order to obtain the leading 
niH contributions. 

The two-loop radiative corrections have the same sign as the one-loop 
ones, and thus result in an enhancement of the Higgs width. The two-loop 
corrections become as large as the one-loop ones for run = 930 GeV. For 
this value of the Higgs mass, the sum of one-loop and two-loop radiative 
corrections is as large as 26% of the tree level widths. This is indicative for 
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the point beyond which the perturbative series stops converging at all in this 
renormalization scheme, and calculations performed by means of Feynman 
diagrams become unreliable. 

This result is rather surprising, considering that for ran = 930 GeV 
the one-loop corrections are quite small, at 13% level. They only become 



substantial for a Higgs as heavy as 1.3 TeV [12|. Considering that most of 
the existing calculations in the electroweak theory were done in the OMS 
renormalization scheme, this raises the question of the validity range of the 
calculations of other processes involving the spontaneous symmetry breaking 
sector, such as the WW — > WW scattering which is of interest in view of 
searches for the Higgs boson at future hadron colliders. 
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Appendix 

Here we give some details related to the use of the techniques of ref. to 
evaluate the scalar integrals needed for this calculation. 

The following two basic functions were introduced in ref. [30]: 



g(mi,m 2 ,m 3 ;k z ) = I dx[Sp(- ) + Sp(- 



f(mi,m2,m 3 ;k 2 ) = I dx 



o 1-2/2 

+2/1 log Vl --. + 2/2 log V2 n } , (15) 
2/1-1 2/2-1 

1-fi 2 



o 2k 

1 2 , 2/1 1 2 , 2/2 i 

-2/1 log 7-7j2/2log T ], 16 

2 2/1 - 1 2 2/2 -1 



where 



1 + K 2 - I? ± yfK 

Vl > 2 = 2^ 

A = (1 + k 2 - ii 2 f + 4kV - 4iK 2 ?/ , (17) 

and 
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No attempt was done to further integrate these functions analytically, 
since this cannot be done in terms of known and easy to calculate functions, 
such as poly logarithms. It is presumably possible to relate them to the Lau- 



ricella functions [17], but it is not clear that this would lead to an efficient 
way to calculate them. Instead, a FORTRAN routine was written to inte- 
grate numerically the g and / functions, as well as the necessary derivatives 
of g, to the desired accuracy. This can be done easily by using an adaptative 
deterministic integration algorithm. 

Some tricks were used to perform the integration in an efficient way. 
First, one extracts the singularities at the ends of the integration path, 
which are of logarithmic type, by a convenient change of variables, such as 
t = \fx. Then, the program chooses the appropriate integration path in the 
complex plane of the integration variable x, in order to avoid the eventual 
singularities of the integrand. The aim is twofold. First, on such a path 
the integrand has a small variance and is smooth, so the integral can be 
calculated to high accuracy by using a small number of points. Second, 
by carrying out the integration on such a path, one avoids automatically 
the numerical instabilities due to large cancellations which occur in the 
computation of the integrand near its branching points. To compute the 
suitable integration path, one needs to know that the integrands of the g 
and / functions given in eq. 15 and 16 have two branching points in the plane 
of the complex Feynman parameter x. They lie on the real axis when the 
functions are calculated above the physical threshold —k 2 > (mi+m2+m^) 2 : 
and their location is given by: 



1 



si,2 = 77~2 [~ a + b + n\ ± J (a - b - fij) 2 - 4fyif] 



Hl t2 = 1-k 2 =f2 v /Z k2 . (19) 
The functions g and / are the finite parts of the following scalar integrals: 
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g(m 1 ,m 2 ,m 3 ;k 2 ) 

1 



J d n pd n q- 



(p 2 + m 2 ) 2 [(q + k) 2 + m 2 ] [(p + c/) 2 + 777.3] 
2 l r . . „ „, , . 1 7T 2 

e 2 ' r _ "" ov i/J ' 4 ' 12 



vr 4 { 3 + - [-1 + 2 7 + 2 log(7T 777?)] + 1 + 

+i[-l + 2 7 + 21og(vr 777 2 )] 2 - 1+^(777!, 777 2 , 7773; A; 2 )} , (20) 

JP"(777i, 777 2 , 7773; & 2 ) = 

_ / dn pd n ^ + ^ = 

J " ^" « (p2 + m 2)[( 9 + fc )2 + m 2 ](r 2 + m 2 )2 

^V{-^ + ^- i[ 7 + log(™ 2 )]+ /(777 1; 7772 , 7773 5 /C 2 )} , (21) 

Let us introduce the following notation: 



G(mi,a 1 ;m2,a2;m 3 ,a 3 ;k 2 ) = 

[ d n vd n a - (221 

J (p 2 + m\ ) ai (q 2 + 777^)«2 [(r + k) 2 + m 2 .]^ ' v ; 

All G scalar integrals can be obtained from Q and T with the help of 
the following relations: 



G(777i,ai + 1; 777-2,02; "73,03; = 
1 Ft 

■G(m 1 ,ai;m2,a2;m3,a 5 ;k ) 



a\ dm 2 
G(t77i,1; 7772,1; 7773, l;k 2 ) = 

-{"7 2 ^(7771, 7772, 777 3 ; k 2 ) + 777-2 ^( m 2 , rU\ , 777 3 ; A: 2 ) 



3 — 77 

+777^ ^(7773, mi, 777 2 ;A: 2 ) +^"(7771, 7772, 7773; A; 2 )} , (23) 

We further notice that any two-loop scalar integral is either a G integral, 
or it can be written as an integral of a certain G after combining all propaga- 
tors with the same loop momentum (p, q or p + q) by introducing Feynman 
parameters. Where further numerical integrations are needed, the numeri- 
cal techniques are similar to those used for the computation of the g and / 
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functions. The main trick is to perform the integration over the remaining 
Feynman parameters on a complex path to avoid the eventual singularities 
of the integrand. The speed of the integration increases dramatically if one 
uses an optimized integration path, along which the integrand is smooth 
enough. The path was defined by means of spline functions. 

Some scalar integrals need to be evaluated at vanishing external mo- 
mentum. In this case T vanishes, and the function g can be integrated 
analytically in terms of Spence functions [14, 3C]: 



g(mi, 1712, m3; 0) = 

1 - i log a log 6- ?L±±Z± [Sp (-^) + Sp(-^) 

1 9 UO 1 , 9 VO 1 9 Ml 1 , n UO IT 2 n , 

+ log 2 ^ + -log 2 ^ + -log 2 ^--log 2 ^ + — , 24 
4 vi 4 ui A vi 4 V2 6 

where 

«i, 2 = ^(1 + 6-aiVA 7 ) 

vi,2 = ~(l-6 + a± V^A 7 ) 
A' = l-2(a + 6) + (a-6) 2 . (25) 

Only two masses appear in this calculation: run, the mass of the Higgs 
boson, and \/f mvy, the mass of the Goldstone modes, with the Landau 
gauge limit taken in the end of the calculation. Therefore one needs some 
mass expansions of the g(mi, rri2, ms; 0) function. 

The necessary expansions are given in the following, sometimes with 
unnecessary precision: 



j(M, M, m) 

j(M, 771, 777 ) 

j(m, M, m) 



1 + k^ ) + x2 -5 + 31ogx +0(x 3, 



3(3 



7T 



+ x(-2 + 21ogx) + x 1 



3 vr 2 

-2 + T + 



+3 log x + log 2 x) + 0(x 3 



(26) 



(27) 



7T" 

1T 



log 2 X 



7T- 



+ x(2 - ^ - 21ogx - log 2 x) + 



16 



+x 2 (y - vr 2 - 71ogx- 31og 2 x) + 0(x 3 ) (28) 

• / , , ,,n „ log 2 x 13 — 6 log x 
j(m,M,M) = -2 + logx-^— +x + 

9 26 — 15 log x , o . , . 

+X 300 + ° {X) (29) 



J [1A01 (M,M,m) = ^[l + x(l-^) + 0(x 2 )] (30) 
jM-O] (M, m,m) = ^2[-2xlogx + C(x 2 )] (31) 

i [1 '°' 0l (m,M,m) = -^[-^^ - 21ogx - -x(tt 2 + 121ogx + 

+31og 2 x) + 0(x 2 )] (32) 
Jinni/ , , 1 r l — logx 7— 6 log x 

+l 37 "30 o'° gX +0 ^l (33) 
jW-»)(M.A/,m) = p[-lt. 2 '^% 0(,')] (34) 

— 7 + 61ogx . 9n1 , N 

+x _^ + 0( a ; 2 )] (35) 

j [0 ' lfi] (M,m,m) = -^2[logx + x(y +41ogx + 

+ log 2 x) + 0(x 2 )] (36) 
j[°' 1 '°](m,m,M) = p[-y -21ogx-log 2 x + x(4-4y - 

-121ogx-41og 2 x) + C(x 2 )] (37) 
j [QXQ \m,M,m) = ^[logx + x(y +71ogx + 

+ log 2 x) + 0(x 2 )] (38) 
j^(m,M,M) = ^Z^l^zL^l^)] (39) 

2 

Here, x = jjp, and the following notations were introduced: 

j(m 1 ,m 2 ,m 3 ) = g(m 1 ,m 2 ,m 3 ;0) - 1 , 
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.[afl 7 ], n Qa & 91 ■( \ 

,1 ^j (roi> ^ > ^ ) = ___ 7 _^_ J ( roi>ro2>ro3) , 

and similar for the derivatives of /. 

One also needs the following relations: 

jW\{M,M,M) = ^-MCl(-)-2j^°^(M,M,M) (40) 

o o 

j(M, M, M) = --J=Ci(|) (41) 

f(M,M,M;-M 2 ) = ~ (42) 

f(M,m,m;-M 2 ) = ~+0(x) (43) 
/[0A0,i] (M;m)m; _ M 2 ) = _l_L + 0(x) (44) 
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Figure captions 



Fig.l The topologies of the one-loop selfenergy diagrams. 

Fig. 2 The topologies of the two-loop selfenergy diagrams. 

Fig. 3 The topologies of the one-loop proper vertex diagrams. 

Fig. 4 The topologies of the two-loop proper vertex diagrams. 

Fig. 5 Triangular diagrams which display endpoint singularities. The 
solid line denotes the Higgs, and the dashed one the Goldstone bosons. 

Fig. 6 The radiative corrections to the partial decay width of the Higgs 
to vector bosons in the one-loop (solid line) and the two-loop (dashed line) 
approximations as a function of the mass of the Higgs boson. 
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